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A MEASUREMENT OF WEAK LENSING BY LARGE SCALE STRUCTURE IN RCS FIELDS 
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ABSTRACT 

We have analysed ~ 24 square degrees of Rc-b&nd imaging data from the Red-Sequence Cluster 
Survey (RCS) , and measured the excess correlations between galaxy ellipticities on scales ranging from 1 
to 30 arcminutes. We have used data from two different telescopes: ~ 16.4 square degrees of CFHT data 
and ~ 7.6 square degrees of CTIO 4-meter data, distributed over 13 widely separated patches. For the 
first time, a direct comparison can be made of the lensing signal measured using different instruments, 
which provides an important test of the weak lensing analysis itself. The measurements obtained from 
the two telescopes agree well. For the lensing analysis we use galaxies down to a limiting magnitude of 
Rc = 24, for which the redshift distribution is known relatively well. This allows us to constrain some 
cosmological parameters. For the currently favored ACDM model (fi m = 0.3, SIa = 0.7, T = 0.21) 
wo obtain <j 8 = 0.811^9 (95% confidence), in agreement with the results from Van Waerbeke et al. 
(2001) which used fainter galaxies (and consequently higher redshift galaxies). The good agreement 
between these two very different weak lensing studies demonstrates that weak lensing is a useful tool in 
observational cosmology. 

Subject headings: cosmology: observations — dark matter — gravitational lensing 



1. INTRODUCTION 

Weak gravitational lensing has proven to be a powerful 
tool for studies of the mass distribution in rich clusters of 
galaxies (for a review see Mellier 1999). Since the pioneer- 
ing work by Tyson, Wenk, & Valdes (1990) much progress 
has been made, and nowadays the weak lensing signals in- 
duced by clusters of galaxies at intermediate redshifts can 
be measured without much difficulty. 

The development of advanced techniques to correct 
for the various observational distortions, such as the 
anisotropy of the point spread function (PSF), the circu- 
larization by the PSF, and the camera induced distortion, 
has been a crucial step, resulting in well calibrated sig- 
nals (e.g., Kaiser, Squires, & Broadhurst 1995; Luppino & 
Kaiser 1997; Hoekstra et al. 1998; Kuijken 1999; Refregier 
2001). Another important development in recent years is 
the advent of mosaic CCD cameras, which enable us to 
quickly image large portions of the sky. 

These advances have made it possible to pursue one of 
the most difficult measurements in the field of weak lens- 
ing: the measurement of the coherent distortions of the 
images of faint galaxies caused by lensing by intervening 
large scale structure, the so-called 'cosmic shear'. The 
analysis of this lensing signal provides an important direct 
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measurement of the statistical properties of the large scale 
mass distribution (e.g., Blandford et al 1991; Kaiser 1992; 
Bcrnardeau, van Waerbeke, & Mellier 1997; Schneider et 
al. 1998). 

Compared to many other methods that are used to con- 
strain cosmological parameters, weak lensing has the ad- 
vantage that no assumptions about the light distribution 
are required. However, weak lensing in itself cannot con- 
strain all the parameters, because of degeneracies between 
them. Better constraints can be obtained by comparing 
weak lensing studies that probe different redshifts (e.g., 
Hu 1999) or when these data are combined with measure- 
ments of the fluctuations of the cosmic microwave back- 
ground (e.g., Hu & Tegmark 1999). 

By now, several groups have reported the detection of 
an excess correlation between galaxy ellipticities, and have 
argued that this signal is caused by lensing by large scale 
structure (e.g., Bacon et al. 2000; Kaiser, Wilson, & Lup- 
pino 2000; van Waerbeke et al. 2000; Wittman et al. 2000; 
Maoli et al. 2001; van Waerbeke et al. 2001). Maoli et 
al. (2001) combined their own results with published mea- 
surements in an attempt to obtain constraints on os, the 
normalisation of the power spectrum, and f2 m , the matter 
density of the universe. They found good agreement with 
studies of cluster abundances. However, the data set stud- 
ied by Maoli et al. (2001) is very inhomogeneous, which 
limits the accuracy of such a direct comparison. 

After these initial detections, which demonstrated the 
feasibility of the method, the obvious next step is to ob- 
tain large uniform data sets. The first results from such 
a survey were presented by van Waerbeke et al. (2001), 
who measured a highly significant lensing signal from 6.5 
square degrees of deep imaging data. 

In this paper we present the results from our analysis of 
~ 24 deg 2 of _Rc7-band data from the Red-Sequence Clus- 
ter Survey (RCS) (e.g., Gladders & Yee 2000), which is a 
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100 deg 2 galaxy cluster survey designed to provide a large 
sample of optically selected clusters of galaxies with red- 
shifts 0.1 < z < 1.4. The data are also useful for a range of 
lensing studies. For example, Gladders, Yee & Ellingson 
(2001) presented the first results for one of the strong lens- 
ing clusters discovered in the survey, for which follow-up 
observations are underway. 

The weak lensing applications are numerous. The sur- 
vey imaging data are relatively shallow compared to what 
is common in weak lensing studies, and as a result the 
statistical uncertainty in the measurements of individual 
structures (such as clusters of galaxies) is large. However, 
thanks to the large survey area, many such structures can 
be detected, and by stacking the signals, one can study 
their ensemble averaged mass distribution (e.g., Hoekstra 
et al. 2001). In addition, follow-up observations will pro- 
vide detailed information on individual systems. 

In this paper we concentrate on the measurement of the 
weak lensing signal induced by large scale structure (cos- 
mic shear). A study of galaxy biasing, based on some 
of these data, is presented in Hoekstra, Yee, & Gladders 
(2001), and a study of the properties of galaxy halos is 
currently underway. 

Compared to other cosmic shear studies, the RCS data 
are shallow, and consequently the signal at a given scale 
is much lower, as is the signal-to-noise ratio. However, 
measuring the weak lensing signal from a shallow survey 
also has several advantages. Down to a limiting magni- 
tude of Rq ~ 24 star-galaxy separation works well. In 
deeper surveys many sources have sizes comparable to the 
size of the PSF, and applying size cuts may change the 
rcdshift distribution of the sources in a systematic way. In 
addition, down to Rc ~ 24 the redshift distribution of the 
sources is fairly well determined. In order to relate the 
observed cosmic shear signal to cosmological parameters, 
a good understanding of the source redshift distribution is 
crucial. 

One worry is the effect of intrinsic alignments of the 
source galaxies, which introduces an additional signal (e.g., 
Heavens et al. 2000; Catelan et al 2001; Crittenden et al. 
2001; Mackey et al. 2001). The amplitude of the effect is 
not well determined, but it is clear that it becomes more 
important for shallower surveys. However, the predictions 
indicate that for a median redshift of z — 0.5 (which is 
similar to our sample of source galaxies) the signal caused 
by intrinsic alignments is still small compared to the lens- 
ing signal (e.g., Makey et al. 2001), and we will ignore the 
effect in this paper. 

In §2 we describe the Red-Sequence Cluster Survey 
(RCS) from which we have used the i?c-band data for 
the analysis presented here. §3 deals with the analysis of 
the data, as well as the corrections for systematic distor- 
tions, such as PSF anisotropy and the distortion by the 
camera. In §4 we discuss the expected signal from weak 
lensing by large scale structure. The results of the analysis 
are presented in §5. 

2. DATA 

2.1. The Red-Sequence Cluster Survey 

The Red-Sequence Cluster Survey 1 (RCS) is a galaxy 
cluster survey designed to provide a large sample of op- 

1 http : //www . astro . utoront o . ca/ ~gladder s/RCS 



tically selected clusters of galaxies with redshifts 0.1 < 
z < 1.4. The planned survey will cover 100 square de- 
grees in both Rc and z' , and consists of 22 widely sepa- 
rated patches of ~ 2.1 x 2.3 degrees. The northern half of 
the survey is observed using the CFH12k camera on the 
CFHT, and the data from the southern half are obtained 
using the Mosaic II camera on the CTIO 4m telescope. 
The patches are imaged down to a 5a point source depth 
of 25.2 magnitudes in the i?c-band, and 23.6 magnitudes 
in the z' filter. 

For the weak lensing analysis we use only the Rc band 
data, as these provide a sufficiently high number density of 
sources to warrant an accurate measurement of the lens- 
ing signal. We present the results based on ~ 16.4 deg 2 of 
CFHT data and - 7.6 deg 2 of CTIO data. In this paper 
we use a subset of the RCS and the data for each patch are 
not contiguous. Thus the largest scale we consider here is 
that of one pointing of the CFH12k or Mosaic II camera. 

We have used data from all 10 patches that have been 
observed using the CFHT, resulting in a total of 53 point- 
ings. The integration times are 900s per pointing. In addi- 
tion we have used three patches (resulting in an additional 
23 pointings) from the first CTIO run, which have integra- 
tion times of 1200s. Some details about the observations 
are listed in Table 1. 

2.2. Data reduction 

Given the large amount of data collected in the sur- 
vey, special attention was paid during the survey design 
on how to handle the data flow. To simplify the construc- 
tion of the science images the data were acquired without 
dithering. Although the gaps between the chips, cosmetic 
defects, and cosmic rays result in a minor loss in area, the 
advantage in handling the data flow is tremendous. The 
loss of area does not affect the result of the weak lens- 
ing analysis, and cosmic rays are easily removed from the 
galaxy catalogs. 

The individual chips from the mosaic imagers are de- 
biased, and flat-fielded using standard techniques. The 
images are used for the object analysis, which is described 
below. A detailed discussion of the reduction pipeline is 
presented in Gladders & Yee (in preparation). 

3. OBJECT ANALYSIS 

Our weak lensing analysis technique is based on that 
developed by Kaiser, Squires & Broadhurst (1995) and 
Luppino & Kaiser (1997), with a number of modifications 
which are described in Hoekstra et al. (1998) and Hoekstra 
et al. (2000). 

This correction scheme assumes that one can model the 
PSF as an isotropic function, convolved with a compact, 
anisotropic kernel. The method does not make any as- 
sumptions about the profile of the PSF or the galaxy, as 
these parameters are measured from the actual data. 

In real data, the PSF is likely to be more complex, and 
the assumption stated above is not valid. However, as 
shown by Hoekstra et al. (1998) matching the measure- 
ment of the PSF parameters to the size of the galaxy re- 
sults in accurate corrections for PSFs with varying ellip- 
ticity as a function of radius (as is the case for the WFPC2 
PSF). 
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Table 1 

(a) Some information for the 53 CFHT pointings used in this 
analysis. Ten widely separated patches have been observed. 
Typically 5 pointings per patch were used, except for the 1417 (4 
pointings) and the 0223 patch (8 pointings); (b) Information for 
the 23 CTIO pointings. We used all the data obtained during 
run 1, and as a result the number of pointings per patch varies. 
The seeing was determined for both telescopes using the sizes of 
stars on chip 3. 



orientations with respect to the PSF. 

We analyse the chips of each pointing separately. After 
the catalogs have been corrected for the various obser- 
vational effects, they are combined into a master catalog 
which covers the observed field (for each pointing). 

The first step in the analysis is to detect the faint galaxy 
images, for which we used the hierachical peak finding al- 
gorithm from Kaiser et al. (1995). We select objects which 
are detected with a significance greater than 5cr over the 
local sky. 

We use single exposures for our analysis, and conse- 
quently cosmic rays have not been removed. However, 
cosmic rays are readily eliminated from the photometric 
catalogs: small, but very significant objects are likely to 
be cosmic rays, or artefacts from the CCD. The peak finder 
gives fair estimates of the object size, and we remove all 
objects smaller than the size of the PSF. 

The objects in this cleaned catalog are then analysed, 
which yield estimates for the size, apparent magnitude, 
and shape parameters (polarization and polarizabilities). 
The objects in this catalog are inspected by eye, in order 
to remove spurious detections. These objects have to be 
removed because their shape measurements are affected 
by cosmetic defects (such as dead columns, bleeding stars, 
halos, diffraction spikes) or because the objects are likely 
to be part of a resolved galaxy (e.g., HII regions). The vi- 
sual inspection is important as it is not possible to remove 
all spurious detections in a fully automatic process. Their 
removal is crucial for an accurate measurement of cosmic 
shear, because they increase the measurement of the vari- 
ance, and introduce artificial ellipticity correlations. 
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In addition, the accuracy of this method has been stud- 
ied in great detail (e.g., Erben et al. 2001; Bacon ct al. 
2001) and the results demonstrate that it works well down 
to the required accuracy for current data sets. Hence, the 
correction scheme is accurate even if the PSF does not 
satisfy the assumption made in the derivation. However, 
this can be understood easily, because any residuals in- 
duced by higher order moments of the PSF are supressed 
as one averages the shapes of galaxies which have random 



Fig. 1. — (a) Plot of the apparent Rp-band magnitude versus 
the half-light radius. The vertical sequence of points at r/j ~ 0'.'45 
(indicated by the rectangle) consists of bright, unsaturated stars, 
that are used to study the PSF. (b) Histogram of the number of 
objects of given ry t . (c) The dashed line shows the fraction of ob- 
jects with half light radii larger than the PSF. The solid line shows 
the fraction of objects larger than the PSF, when the total counts 
are corrected for the contribution by stars. Objects larger than the 
PSF arc assumed to be galaxies and only these are used in the weak 
lensing analysis. The figure demonstrates that down to R ~ 24 
this separation selects most galaxies, as almost all objects arc larger 
than the PSF. 
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3.1. Correction for the PSF 

To measure the small, lensing induced distortions in the 
images of the faint galaxies it is important to accurately 
correct the shapes for observational effects, such as PSF 
anisotropy and seeing; PSF anisotropy can mimic a cosmic 
shear signal, and a correction for the seeing is required to 
relate the measured shapes to the real lensing signal. 

To do so, we follow the procedure outlined in Hockstra 
et al. (1998). We select a sample of moderately bright 
stars from our observations, and use these to characterize 
the PSF anisotropy and seeing. Figure la shows a plot of 
the apparent i?c-band magnitudes of the detected objects 
versus their measured half-light radii for one of the chips 
of the A2 pointing of the 0223 patch (seeing ~ Of! 77). We 
have also indicated the region from which we have taken 
the stars used for the analysis of the PSF. 

We fit a second order polynomial to the shape parame- 
ters of the selected stars for each chip. This procedure is 
repeated for various dispersions of the weight function (for 
details see Hoekstra et al. 1998). In the left panel of Fig- 
ure 2 we present the resulting model PSF anisotropy for 
the A2 pointing of the 0223 patch. To show in more detail 
the higher order spatial dependence of the anisotropy we 
have subtracted the average ellipticity. Although the fits 
were obtained from the individual chips, the mosaic image 
in Figure 2 shows continuity between the chips. 

The results for one of the CTIO pointings is presented 
in the right panel of Figure 2. Comparison of the patterns 
presented in Figure 2 with other pointings shows that the 
pattern is fairly stable, although the amplitute varies, be- 
cause of focus variations. In general the PSF anisotropy is 
small, a point which we will address in more detail below, 



when we examine the residuals left after correction of the 
shapes of the galaxies used in the weak lensing analysis. 

3.2. Telescope distortion 

The effect of the PSF is not the only observational dis- 
tortion that has to be corrected. The optics of the camera 
stretches the images of galaxies (i.e., it introduces a shear) 
because of the non-linear remapping of the sky onto the 
CCD. We have used observations of astrometric fields to 
find the mapping between the sky and the CCD pixel co- 
ordinate system, and derived the corresponding camera 
shear. 

The camera shear for the CFH12k camera for run 2 is 
presented in Figure 3. The shear introduced by the cam- 
era is small, reaching a maximum value of ~ 1% at the 
edges of the field of view. 

Other weak lensing studies, which use dithered observa- 
tions, need to remap the images before combining the data, 
thus removing the camera distortion. We have analysed 
single exposures, and measured the shapes of the galaxies 
on the reduced images that have not been remapped to 
remove the camera distortion. As discussed in Hockstra 
et al. (1998), the images of both the stars and galaxies are 
sheared by the camera. The measured shape of the PSF 
(as shown in Figure 2) is then the combination of PSF 
anisotropy and camera shear, and therefore the real PSF 
anisotropy is somewhat smaller than the measured value. 
The change in the ellipticity of an object caused by the 
camera shear depends on its shape (and hence the cor- 
rection) whereas the PSF anisotropy correction depends 
mainly on the size of the object. As a result the correction 
for PSF anisotropy leaves a residual ellipticity, because the 
correction will be too large for larger galaxies. However, 
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Fig. 2. — Left panel: PSF anisotropy as a fuction of position for one of the CFHT pointings; right panel: PSF anisotropy for one of the 
CTIO pointings. One of the chips in the CTIO observations was not functioning and has been omitted. The sticks indicate the direction of 
the major axis of the PSF, and the length is proportional to the observed ellipticity of the PSF. In order to show the higher order spatial 
dependence of the anisotropy we have subtracted the average ellipticity. The direction of the average PSF anisotropy is indicated in the top 
right box, and the amplitude is indicated as well. Although the PSF anisotropy was determined from fits to the observed shapes for individual 
chips, the figure clearly shows a large scale dependence. 
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Fig. 3. — Left panel: camera distortion measured for the CFH12k camera for run 2. The shear introduced by the camera is small, reaching 
a maximum value of ~ 1% at the edges of the field of view; Right panel: the camera distortion measured from the CTIO data. One of the 
chips was not functioning and has been omitted. The induced shear by this camera is somewhat larger than the CFH12k camera, but still 
small: ~ 2% at the edges of the field of view. 



Hockstra et al. (1998) demonstrated that the correction 
for the residual camera shear is straightforward: one just 
needs to subtract the camera shear from the measured 
shear, which is what we have done. 

The camera shear is more or less radial with respect 
to the center of the camera (although the camera shear 
of the CTIO Mosaic II camera shows a significant non- 
radial component), which results in a negative tangential 
shear. It is therefore useful to examine the average tan- 
gential shear of the galaxies with respect to the center of 
the camera. The results are presented in Figure 4a. The 
solid circles indicate the average tangential distortion of 
the galaxies with respect to the center of the CFH12k cam- 
era after correction for PSF anisotropy. These measure- 
ments agree well with inferred camera shear (solid line). 
After subtracting the camera shear we obtain the open 
circles, which are consistent with no signal. 

3.3. Residuals 

The correction scheme developed by Kaiser et al. (1995) 
has been tested extensively (e.g., Hoekstra et al. 1998; Ba- 
con et al. 2001; Erben et al. 2001). The assumptions that 
have to made to derive the original Kaiser et al. (1995) cor- 
rection parameters do not necessarily hold in real data, the 
modifications suggested by Hoekstra ct al. (1998) allow it 
to be applied to more complicated PSFs. This is supported 
by numerous simulations which indicate the method works 
remarkably well down to the required accuracy for current 
data sets. 

In addition we have tested the method using a realistic 
simulation. The simulated data sets were created using 
the software tools SkyMaker and Stuff 2 , which have been 
described in detail in Erben et al. (2000). The simulated 
galaxies have realistic profiles, with a mix of early type, 
late type galaxies, and disk/bulge ratios matched to actual 
observations. The PSF is computed using realistic pupil 

2 see http://terapix.iap.fr/soft 



functions, and includes all the problems encountered in 
real data, such as coma, tracking errors, aberration, spi- 
der arms from the support of the secondary mirror, etc. 

The simulation is described by Van Waerbeke et al. (in 
preparation), and here we briefly discuss the results. Van 
Waerbeke used an N-body simulation to infer the corre- 
sponding (cosmic) shear by means of ray-tracing. 
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Fig. 4. — (a) Solid circles indicate the average tangential distor- 
tion of the galaxies with respect to the center of the CFH12k camera 
after correction for PSF anisotropy. The solid line corresponds to 
the average tangential camera shear. The open circles give the mea- 
surements of the galaxies after correcting for both PSF anisotropy 
and camera distortion. The results indicate that the two steps in the 
correction have worked well, (b) number of galaxies as a function 
of radius used to produce panel (a). 



6 



Weak lensing by large scale structure 




e, (star) e a (star) 

Fig. 5. — (a) (71) for the galaxies with 20 < R < 24 as a function of the value of the ellipticity component ei of the stars used to correct 
the galaxies. The results show the expected strong correlation; (b) The average 71 after correction for PSF anisotropy (note the vertical scale 
has been expanded with respect to panel a); (c) The number of galaxies with a given value for ei of the PSF; (d)-(f) same as (a)-(c), but for 
72. The residuals are consistent with no signal, demonstrating that the correction for the PSF anisotropy has worked well. 



Artifical images of galaxies were sheared using these re- 
sults. These galaxies were "observed" and the final images 
are convolved with a realistic PSF. A mosaic of 900 im- 
ages (corresponding to <~ 12 deg 2 ), each with a different 
PSF was given to us. We note that the PSFs used in the 
simulation had worse systematics than the once observed 
in the RCS data. 

We analysed these images in the same way as we have 
analysed the RCS data, and measured the cosmic shear sig- 
nal. In doing so, we were able to recover the input lensing 
signal, which was unknown to us. We were able to recover 
the signal within 10% (<~ la) of the input value. Because 
of the noise introduced by the intrinsic shapes of the galax- 
ies, larger simulations are required to test whether we can 
measure the lensing signal to even higher accuracies. 

We have also examined the residuals in the RCS data af- 
ter PSF correction. A useful test, although not definitive, 
is to plot the average shape of the galaxies as a function of 
the shape of the PSF. The results for the two components 
of the shear are presented in Figure 5. Panels (c), and (f) 
indicate that the PSF anisotropy is small for most galax- 
ies. The residuals presented in panels (b), and (e) (note 
the different vertical scale) are consistent with no signal. 
The results of Figure 5 and the results from the simulation 
suggest that we are able to measure the lensing signal to 
an accuracy better than 10%. 

4. COSMIC SHEAR SIGNAL 

To study the weak lensing signal caused by large scale 
structure, we use the top-hat smoothed variance of the 
shear (Bacon et al. 2000; Kaiser et al. 2000; Maoli et 



al. 2001; van Waerbeke et al. 2000; van Waerbeke et al. 
2001). Other statistics, such as the ellipticity correlation 
function (Kaiser 1992; van Waerbeke 2001; Wittman et 
al. 2000) or the aperture mass statistic (Schneider 1998, 
Schneider et al. 1998; van Waerbeke 2001) have also been 
used. 

Here we briefly discuss how the lensing signal depends 
on the assumed cosmology and the redshift distribution of 
the sources. Detailed discussions on this subject can be 
found elsewhere (e.g., Schneider et al. 1998; Bartelmann 
& Schneider 2001). 

Given a cosmological model, the variance in the shear 
caused by large scale structure can be computed as a func- 
tion of aperture size 9 (e.g., Jain & Seljak 1997): 



<7 2 >W 



= 2tt ( 
Jo 



dl I P K {1) 



Ji(W) 



ttW 



(1) 



where 9 is the radius of the aperture used to compute the 
variance, and Ji is the first Bessel function of the first 
kind. The effective convergence power spectrum P K (l) is 
given by 



Pk{1) 



4c 4 



I-Wh / 

I n 



W(w) 
a(w) 



I 



fK{w)' 



(2) 

Here w is the radial coordinate, wh corresponds to the 
horizon, a(w) the cosmic scale factor, and Jk{w) the co- 
moving angular diameter distance. As shown by Jain & 
Seljak (1997) and Schneider et al. (1998) it is necessary to 
use the non-linear power spectrum in equation (2). This 
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power spectrum is derived from the linear power spectrum 
following the prescriptions from Peacock & Dodds (1996). 

W(w) is the source-averaged ratio of angular diame- 
ter distances Di s /D s for a redshift distribution of sources 
Pw(w): 



W(w) 



a i i '\JK(w - w) 

dwPw{w) feK) ■ (3) 



Thus it is important to know the redshift distribution 
of the sources, in order to relate the observed signal to 
P K (l). A detailed discussion of the adopted redshift dis- 
tribution can be found in Section 4.1. Figure 6 shows 
the top- hat smoothed variance (j 2 ) on a scale of 1 ar- 
cminute as a function of limiting magnitude of the sample 
of sources. To compute the signal we used galaxies with 
20 < R < -Rii m and used the photometric redshift distribu- 
tion inferred from the Hubble Deep Fields north and south 
(see section 4.1). The top axis indicates the corresponding 
median redshift of the source galaxies. 

One of the advantages of deep observations is obvious 
from Figure 6: the signal increases quickly for limiting 
magnitudes fainter than R = 24 (or higher median red- 
shift). In addition the number of sources increases as well, 
resulting in higher signal-to-noise ratios of the measure- 
ments. 



4.1. Redshift distribution 

In order to relate the observed lensing signal to physical 
parameters, such as as or Q m , knowledge of the redshift 
distribution of the sources is crucial. The galaxies used in 
weak lensing surveys are generally too faint to be included 
in redshift surveys, and little is known about their redshift 
distribution from spectroscopic studies. 

Compared to the other, deeper, cosmic shear studies, 
our analysis has the major advantage that the redshift 
distribution of the sources we use is better known. Down 
to a limiting magnitude of Rc = 24 the redshift distribu- 
tion has been determined spectroscopically by Cohen et al. 
(2000), although this survey is limited to a relatively small 
patch of sky, and is likely to suffer from cosmic variance. 

In addition the galaxies are larger, which is demon- 
strated in Figure 1: down to Rc = 24 the galaxies are 
easily separated from the stars. This has the advantage 
over deeper surveys (where the fainter galaxies have sizes 
comparable to the PSF) in that selecting objects larger 
than the PSF does not change the redshift distribution 
significantly. Even for the worst seeing images considered 
here (~ l'.'l) the stars are well separated from the galaxies 
for R c = 24. 

Photometric redshift studies, in particular those based 
on the Hubble Deep Fields (e.g., Fcrnandez-Soto et al. 
1999; Chen et al. 1998) have also provided important in- 
formation. The results of Hoekstra et al. (2000) have 
demonstrated that they generally work well. However, 
Hoekstra et al. (2000) noted a difference between the red- 
shift distributions inferred for the Northern and the South- 
ern field, and such field to field variation is not unexpected. 
However, currently little is known about the amplitude of 
such variations, and more studies are required to constrain 
the redshift distributions of these faint galaxies. 
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Fig. 6. — The expected variance induced by large scale struc- 
ture in an aperture of radius 1 arcminutc, as a function of lim- 
iting magnitude of the source galaxies. To derive these results 
we used source galaxies fainter than R = 20, and assumed per- 
fect shape measurements. The result for the SCDM (solid line, 
flm = 1, Q A = 0, a$ = 0.5, r = 0.5), OCDM (dashed line, 
n m = 0.3, Q A = 0, as = 0.85, T = 0.21) and ACDM (dotted 
line, U m = 0.3, U A = 0.7, a s = 0.9, T = 0.21) arc shown. The 
cosmic shear signal increases rapidly with increasing limiting mag- 
nitude, or increasing median redshift (as indicated by the top axis). 



In order to minimize the contribution of cosmic vari- 
ance to the redshift distribution, we use the photometric 
redshift distributions from Fernandez-Soto et al. (1999) 
and Chen et al. (1998) to compute the predicted lens- 
ing signal for a given cosmology (see Fig. 2 from Hoekstra 
2001 for the resulting redshift distribution). To do so, 
we have to take into account that the uncertainty in the 
shape measurements depend on the apparent magnitudes 
(and thus on the rcdshifts) of the sources: the contribution 
of distant, small faint galaxies (with noisy shape measure- 
ments) to the measured lensing signal is smaller compared 
to brighter galaxies. 

This is illustrated in Figure 7. Figure 7a shows the ex- 
pected (based on modeling of deep number counts), the 
observed (galaxies for which shapes could be measured), 
and effective number counts (dotted line) as a function of 
apparent magnitude. The effective number density takes 
into account the uncertainty in the galaxy shapes, and 
gives a good indication of which galaxies contribute most 
to the measurement of the lensing signal. Based on the 
results displayed in Figure 7, we decided to use galaxies 
with 20 < Rc < 24 for the lensing analysis. 

The relative weight (normalised to unity for bright 
galaxies) as a function of magnitude is shown in Figure 7b. 
This weight function is simply the inverse square of the 
uncertainty in the shape measurement (see Hoekstra et 
al. 2000 for details), and reflects the fact that the shapes 
estimates of faint galaxies are more noisy. We derive the 
"effective" redshift distribution using this weight function. 
This effective redshift distribution is used to compute the 
predicted lensing signals discussed in Section 5. 
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R [mag] 

Fig. 7. — (a) Observed number density of galaxies (i.e. galaxies 
for which shapes could be determined) as a function of apparent 
Rc magnitude (solid line). The expected number density (based 
on modeling of observed number counts) is indicated by the dashed 
line, (b) The weight (based on the uncertainty in the shape mea- 
surements; see Hoekstra et al. 2000) multiplied by the complete- 
ness fraction as a function of apparent magnitude. The product of 
the number of galaxies and the weight gives a good indication of 
the relative contribution to the lensing signal (dashed line in panel 
(a). The result shows that most of the signal comes from galaxies 
around Rq = 23.5. The vertical dashed lines indicate the region 
20 < Rc < 24, the range of apparent magnitudes for the source 
galaxies we will use in the lensing analysis. 



number of random realisations of the data set where the 
orientations of the galaxies were randomized. 

Note that the measurements at various scales are 
strongly correlated, and this causes all the CTIO mea- 
surements at large scales to be higher than the CFHT 
results. The results obtained from the two different tele- 
scopes agree very well with one another. Because the sys- 
tematics for the two data sets are different, this excellent 
agreement demonstrates that the various observational bi- 
ases have been removed successfully. 

Another useful test is to compare the signals from the in- 
dividual patches. Figure 9 displays the top-hat smoothed 
variance of the shear for an aperture of radius 2.5 arcmin- 
utes (where the signal-to-noise ratio is highest) for the 13 
observed patches. The measurements for the individual 
patches agree with the ensemble average. 

It is also important to examine the correction for the 
circularization by the PSF, because it determines the am- 
plitude of the signal. To do so, we computed the variance 
on a scale of 2.5 arcminutes for each pointing, and looked 
for a correlation with seeing. The results show no trend 
with seeing. 

Figure 10 shows the measurement of the top-hat 
smoothed variance for the full weak lensing data set . The 
signal-to-noise ratio of our measurements is very good, 
reaching a maximum of ~ 6 at a radius of 2.5 arcminutes. 

For comparison the predictions for three different cos- 
mological models are also shown in this figure. All three 
models are good fits to the data, indicating the need for 
additional observational constraints (cf. §5.2) 



5. RESULTS 

5.1. Observed signal 

In this section we present the measurement of the weak 
lensing signal caused by large scale structure using the 
top-hat smoothed variance of the shear. This statistic has 
been used by other groups to detect the cosmic shear sig- 
nal (e.g., Bacon et al. 2000; Kaiser et al. 2000; Maoli 
et al. 2001; Van Waerbeke et al. 2000, 2001). The top- 
hat smoothed variance is fairly insensitive to errors in the 
analysis because residual shears are added in quadrature. 
Consequently the observed signal can always be considered 
as an upper limit, because residual errors always increase 
the variance. However, the results presented in Figure 5 
and the simulations discussed in section 3.3 indicate that 
we can measure the shapes of the galaxies accurately. 

As was found by van Waerbeke et al. (2000), close pairs 
of galaxies can introduce an excess signal, because of over- 
lapping isophotes. We therefore remove pairs with a sep- 
aration of less than 2'.'16 (which corresponds to 10 pixels 
for the CFHT data, and 8 pixels for the CTIO data). This 
lowered the signal at small scales (~ 20% for an aperture 
of radius 1 arcminute). We note that on these scales in- 
trinsic alignments can also be important. To compute the 
top-hat smoothed variance we use the practical estimators 
given in Van Waerbeke et al. (2001). 

Figure 8 shows the top-hat smoothed variance as a 
function of scale for both the CFHT data (filled circles) 
and the CTIO data (open circles) using galaxies with 
20 < Rc < 24. The errorbars are estimated from a large 
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Fig. 8. — Top-hat smoothed variance of the shear as a function of 
aperture radius. The filled circles indicate the measurements based 
on 16.4 deg of CFHT data, and the open circles correspond to the 
analysis of 7.6 deg 2 of CTIO data. For display purposes the CTIO 
points have been offset slightly in radius. Note that the measure- 
ments at various scales are strongly correlated, and this causes all 
the CTIO measurements at large scale to be higher than the CFHT 
results. The results obtained from the two different telescopes agree 
well with one another. 
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FlG. 9. — The top-hat smoothed variance of the shear for an 
aperture of radius 2.5 arcminutes for the 13 observed patches. The 
shaded region corresponds to the la limits around the average. The 
measurements for the individual patches agree well with the average, 
indicating that the cosmic variance is small. 

It is difficult to compare our measurements directly to 
most other studies, because of the difference in the source 
rcdshift distribution (caused by the diffence in filters and 
integration times). However, we can compare directly to 
the results from Bacon et al. (2000), who have used a sim- 
ilar cut in apparent magnitude. Bacon et al. (2000) find 
a variance of (j 2 ) — (24 ± 7) x 10~ 5 in 8 x 8 arcminute 
cells. This scale is similar to a radius of ~ 4 arcminutes for 
which we find (j 2 ) = (2.8 ±0.5) x 10 -5 . Our signal is much 
lower than their result. Bacon et al. (2000) note a residual 
correlation between the shape of the PSF and the galax- 
ies (their Fig. 7), and this might explain their increased 
variance. We note that because of the large errorbar on 
the measurement of Bacon et al. (2000) the results are 
consistent at the 3a level. 



5.2. Constraints on cosmological parameters 

As described in Section 4 the predicted amplitude of the 
top-hat smoothed variance depends on the various cos- 
mological parameters, and therefore provides a powerful 
method to constrain these parameters. Unfortunately sev- 
eral degeneracies exist (e.g., Bernardeau et al. 1997; Jain 
& Seljak 1997). These studies show that the amplitude 
of the signal is mainly determined by a combination of as 
and fi m , although the shape parameter T (<~ fl m h m a 
CDM cosmology) is also important. 

As demonstrated by Van Waerbeke et al. (2001), it 
is possible to partially break the degeneracy between a$ 
and Q TO if priors on the shape of the power spectrum are 
assumed. The value of Qa can be constrained by com- 
bining the measurements of the lensing signal for sources 
at different redshifts, using the fact that the angular di- 
ameter distances depend on Qa- The combination of the 
weak lensing measurements and the constraints from stud- 
ies of the fluctuations of the cosmic microwave background 
(CMB) will provide much stronger constraints of the pa- 
rameters, as the combination will break the degeneracies. 
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Fig. 10. — Measurement of the top-hat smoothed variance (excess 
variance caused by lensing by large scale structure) using galaxies 
with 20 < R c < 24. The data consist of 16.4 deg 2 of CFHT data 
and 7.6 deg 2 of CTIO data. The drawn lines correspond to the best 
fit SCDM (Q m = 1, Q A = 0, T = 0.7, a$ = 0.31; solid line), the best 
fit OCDM (fi m = 0.3, n A = 0.0, r = 0.21, a 8 = 0.86; dashed line), 
and the best fit ACDM (Cl m = 0.3, Q A = 0.7, T = 0.21, <r s = 0.81; 
(dotted line) model. We fixed h = 0.7, which gives a high value 
for r for the SCDM model. Without additional constraints on the 
cosmological parameters, the lensing results are consistent with a 
wide range of cosmological models. The errorbars arc estimated 
from a large number of random realisations of the data set where 
the orientations of the galaxies were randomized. Note that the 
points at various scales are strongly correlated. Under the assump- 
tion that the lensing structures are halfway between the observer 
and the sources, we have indicated a scale of 1 Mpc. 

Here we will use the measurements of the top-hat 
smoothed variance to find constraints on f2 m and as- We 
assume that the measurements follow a normal distribu- 
tion. We computed the covariance matrix from a large 
number of random realisations of the data, thus including 
the survey geometry in the noise correlation. 

We compute the model predictions using equation 3, un- 
der the assumption that T — 0.21, and using the effective 
redshift distribution discussed in section 4.1. We use this 
value for T to allow for a direct comparison with the re- 
sults from Van Waerbeke et al. (2001). The predictions 
are compared to the observations, and the likelihood for 
the combination of O m and a s is computed. The results 
for models with Q A — are presented in Figure 11. The 
contours indicate the 68.3%, 95.4%, and 99.7% confidence 
limits on two parameters jointly. 

The results for the Q m + fl\ = 1 model are presented 
in Figure 12. The best fit value for Q TO is lower than for 
the OCDM model, but no strong constraints on Cl m and 
as can be placed. Allowing for larger values (r ~ 0.7), the 
goodness of fit for high Cl m models is comparable to those 
of low Q TO models. We note, however, that studies of the 
galaxy correlation function suggest values for T = 0.1 — 0.3 
(e.g., Dodelson et al. 2001) 

The best fit ACDM cosmology yields a 8 = 0.81^ ^ 
(95% confidence; T = 0.21). For the best fit OCDM model 
we obtain a s = Q.86~^q']% (95% confidence). Van Waerbeke 
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et al. (2001) find a 8 = 0.99±g;?g (95% confidence) for an 
OCDM model, which is in fair agreement with our result, 
in particular when the uncertainty in the redshift distri- 
bution used by Van Waerbeke et al. (2001) is taken into 
account. 

To investigate the agreement between the CFHT and 
CTIO data we also computed the best fit values for as for 
the two data sets. We obtain a 8 = 0.84 ± 0.09 (68% con- 
fidence) for the CFHT data, and a 8 = 0.76 ± 0.12 (68% 
confidence) for the CTIO data. The values agree well, and 
the probability of finding a smaller difference is 40%. 



0.5 




Fig. 1 1 . — Likelihood contours as a function of Q m and as , inferred 
from the analysis of the top-hat smoothed variance. We have used 
only the measurements at radii > 1 arcminute. The contours have 
been computed by comparing the measurements to CDM models 
with n = 1, T = 0.21 and S7a = 0. The contours indicate the 68.3%, 
95.4%, and 99.7% confidence limits on two parameters jointly. 
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The RCS probes a different redshift range than the 
study presented by Van Waerbeke et al. (2001). Although 
both results are based on the same correction scheme, 
there are many differences in the various steps in the analy- 
ses. Thus the agreement found here suggests that accurate 
measurements of cosmic shear can be made. 

6. CONCLUSIONS 

We have analysed ~ 24 square degrees of i?c-band imag- 
ing data from the Red-Sequence Cluster Survey to study 
the weak lensing signal caused by intervening large scale 
structure. To minimize the effect of cosmic variance, the 
measurements have been obtained from 13 patches that 
are widely separated on the sky. We have used data from 
two different telescopes: ~ 16.4 square degrees of CFHT 
data and <~ 7.6 square degrees of CTIO data. We detect 
the signal with high confidence on scales ranging from 1 
to 30 arcminutes using galaxies with 20 < Rc < 24. 

Because of the various observational distortions which 
need to be corrected for, a careful analysis of the residuals 
is important. The results suggest that we have successfully 
corrected for the systematics. In addition we compared the 
measurements from CFHT and CTIO and find excellent 
agreement. 

Compared to other studies of cosmic shear, the RCS 
imaging data is relatively shallow. This has the disadvan- 
tage that the lensing signal is low. However, the galax- 
ies are larger (which results in smaller corrections for the 
PSF), and the redshift distribution is known fairly well 
(which is important for determining cosmological param- 
eters). 

Intrinsic alignments of galaxies contaminate the lensing 
signal. This is particularly important for lensing studies 
that use low redshift galaxies. The median redshift of our 
sample of sources is <~ 0.5, and predictions of the contri- 
bution of the intrinsic alignments suggest it is small for 
the results presented here. However, the effect of intrinsic 
alignments is expected to be comparable to the errorbars 
of the full RCS data set, and eventually needs to be cor- 
rected for. The RCS will be complemented with B and V 
band imaging which will provide photometric redshifts for 
a large number of galaxies. With such a data set we will 
be able to measure the intrinsic alignments by selecting 
galaxies at similar photometric redshifts. 

We use the photometric redshift distribution inferred 
from the Hubble Deep Fields to relate the measured top- 
hat smoothed variance to estimates of cosmological pa- 
rameters. Because of degeneracies in the parameters we 
can only place constraints on Q m and a 8 jointly. For 
the currently favoured ACDM model (Q TO = 0.3, Oa = 
0.7, T = 0.21) we obtain a s = 0.8lt^g (95% confidence), 
in good agreement with the results from Van Waerbeke et 
al. (2001). 

The RCS data and the observations used by Van Waer- 
beke are quite different, and also the weak lensing analyses 
are somewhat different. Thus the agreement found here 
suggests that accurate measurements of cosmic shear can 
be made. 



Fig. 12. — Same as Figure 12, but now for an Sl m + J1a = 1 CDM 
cosmology. 



It is a pleasure to thank Ludo van Waerbeke for many 
discussions, and providing a useful test of the correction 
scheme used in this paper. We would like to thank the 



Hockstra ct al. 



11 



CFHT and CTIO TACs for their generous allocation of 
telescope time. This research is partially supported via an 
operating grant from NSERC to HKCY. MDG acknowl- 
edges support from NSERC via PGSA and PGSB scholar- 



ships. PBH acknowledges financial support from Chilean 
grant FONDECYT/ 1010981. LFB's research is funded by 
FONDECYT, Chile under Proyecto 1000537. 



REFERENCES 



Bacon, D., Refregier, A., & Ellis, R.S. 2000, MNRAS, 325, 1065 
Bacon, D., Refregier, A., Clowe, D., & Ellis, R.S. 2001, MNRAS, 
325, 1065 

Bartelmann, M., & Schneider, P. 2001, Physics Reports, 340, 291 
Bernardeau, F., van Waerbeke, L., Mellier, Y. 1997, A&A, 322, 1 
Blandford, R.D., Saust, A.B., Brainerd, T.G., & Villumsen, J.V. 

1991, MNRAS, 251, 600 
Catclan, P., Kamionkowski, M., & Blandford, R.D. 2001, MNRAS, 

320, L7 

Chen, H.-W., Fernandez-Soto, A., Lanzetta, K.M., Pascarelle, 

S.M., Puetter, R.C., Yahata, N., & Yahil, A., preprint, 

astro-ph/9812339 
Cohen, J.G., Hogg, D.W., Blandford, R.G., Cowic, L.L., Hu, E., 

Songaila, A., Shopbell, P., & Richberg, K. 2000, ApJ, 538, 29 
Crittenden, R.G., Natarajan, P., Pen, U.-L., & Thcuns, T. 2001, 

ApJ, submitted, astro-ph/0009052 
Dodelson, S. et al. 2001, ApJ, preprint, astro-ph/0107421 
Erbcn, T., Van Waerbeken, L., Bertin, E., Mellier, Y., & Schneider, 

P. 2001, A&A, 368, 766 
Fernandez-Soto, A., Lanzetta, K.M., & Yahil, A. 1999, ApJ, 513, 34 
Gladdcrs, M.D., & Yee, H.K.C. 2000, to appear in "The New Era 

of Wide-Field Astronomy", astro-ph/0011073 
Gladders, M.D., Yee, H.K.C, & Ellingson, E. 2001, AJ, accepted 
Heavens, A., Refregier, A., & Heymans, C. 2000, MNRAS, 319, 649 
Hoekstra, H., Franx, M., Kuijken, K., & Squires, G. 1998, ApJ, 504, 

636 

Hoekstra, H., Franx, M., Kuijken, K. 2000, ApJ, 532, 88 



Hoekstra, H., Franx, M., Kuijken, K., Carlberg, R.G., Yee, H.K.C, 
Lin, H., Morris, S.L., Hall, P.B., Patton, D.R., Sawicki, M., & 
Wirth, CD. 2001a, ApJ, 548, L5 

Hoekstra, H. 2001, A&A, 370, 743 

Hoekstra, H., Yee, H.K.C, & Gladders, M.D. 2001b, ApJ, 558, Lll 

Hu, W. 1999, ApJ, 522, L21 

Hu, W. & Tegmark, M. 1999, ApJ, 514, L65 

Jain, B., & Seljak, U. 1997, ApJ, 484, 560 

Kaiser, N. 1992, ApJ, 388, 272 

Kaiser, N., Squires, C, & Broadhurst, T. 1995, ApJ, 449, 460 
Kaiser, N., Wilson, C, & Luppino, G.A. 2000, ApJL, submitted, 

astro-ph/0003338 
Kuijken, K. 1999, A&A, 352, 355 
Luppino, G.A., & Kaiser, N. 1997, ApJ, 475, 20 
Mackey, J., White, M., & Kamionkowski, M. 2001, MNRAS, 

submitted, astro-ph/0106364 
Maoli, R., Van Waerbeke, L., Mellier, Y., Schneider, P. Jain, B., 

Bernardeau, F., Erben, T. & Fort, B. 2001, A&A, 368, 766 
Mellier, Y. 1999, ARA&A, 37, 127 
Peacock, J.A., & Dodds, S.J. 1996, MNRAS, 780, L19 
Refregier, A. 2001, submitted to MNRAS, astro-ph/0105178 
Schneider, P. 1998, ApJ, 498, 43 

Schneider, P., van Waerbeke, L., Jain, B., & Kruse, G. 1998, 

MNRAS, 296, 873 
Tyson, J.A., Wcnk, R.A., & Valdes, F. 1990, ApJ, 349, LI 
Van Waerbeke, L., et al. 2000, A&A, 358, 30 
Van Waerbeke, L., et al. 2001, A&A, 374, 757 
Wittman, D.M., Tyson, J. A., Kirkman, D., Dell' Antonio, I., & 

Bernstein, G. 2000, Nature, 405, 143 



